Non-adiabatic time-dependent density functional theory of the impurity resistivity of 

metals 
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We make use of the time-dependent density functional theory to derive a new formally exact 
expression for the dc resistivity of metals with impurities. This expression takes fully into account 
the dynamics of electron-electron interactions. Correction to the conventional T-matrix (phase- 
shifts) theory is treated within hydrodynamics of inhomogeneous viscous electron liquid. As a 
first application, we present calculations of the residual resistivity of aluminum as a function of the 
atomic number of the impurities. We show that the inclusion of many-body corrections considerably 
improves the agreement between theory and experiment. 
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Scattering of carriers by impurities is one of the fun- 
damental mechanisms of resistivity in solids [1]. Within 
the single-particle theory, this problem can be efficiently 
addressed with the use of well established techniques of 
potential scattering theory. Electron-electron interac- 
tions can be included, at this level, within the ground- 
state density-functional theory [2J, |3j as electrostatic and 
exchange-correlation effective static potential which scat- 
ters the electrons as single particles |4, 5]. However, the 
single-particle approach fails to account for the dynamical 
exchange and correlation effects which cannot be forced 
into the mold of a static mean field theory. 

A powerful theoretical tool has been devised to ac- 
count, in principle exactly, for dynamical electron- 
electron interaction effects in inhomogeneous systems. 
This is known as the time-dependent density functional 
theory (TDDFT) Q. Historically, TDDFT was devel- 
oped to improve the calculation of atomic scattering 
cross-sections and excitation energies in both bounded [7( 
and extended systems [a |9|. However, in recent years 
this theory has often been applied to the treatment of 
static phe nomena (e.g., the polarizability of polymer 
chains [id Hlj). an d steady-state transport phenomena 
(e.g., the stopping power of metals for slow ions 12 .H3 . 
and the conductance of quantum point contacts Llf- 18j ) . 
In these applications one needs to take the zero-frequency 
limit of a time-dependent process, which can be properly 
described by TDDFT or its current generalization - the 
time-dependent current density functional theory (TD- 
CDFT) [13, Hijj]. 

Moving in the same direction, in this Letter we present 
the first complete TDDFT formulation for the impurity 
resistivity of metals. We derive an exact formula for 
the frequency-dependent resistivity in terms of quantities 
that can be calculated entirely within density functional 
theory. While the standard Kubo formula gives a for- 
mally exact expression for the conductivity, our formula 
gives an expression for the resistivity, and furthermore 
does not require that we calculate explicitly the current 
distribution. A major advantage of working with the re- 



sistivity rather than with the conductivity is that phys- 
ically distinct dissipative processes enter the resistivity 
as additive contributions (Matthiessen's rule [21|). In 
particular, we find that our expression for the resistiv- 
ity naturally separates into a single-particle contribution 
and a dynamical many-body contribution. The former is 
shown to reduce to the classical potential-scattering for- 
mula for the resistivity; the latter takes into account the 
viscosity of the electron liquid 221 ] . 

Our formula appears to be a promising tool for a sys- 
tematic improvement on the existing calculations of the 
resistivity of metals. We demonstrate this in a concrete 
application, namely the model calculation of the resistiv- 
ity of aluminum in the presence of random impurities. 
We show that dynamical corrections, calculated with the 
help of the available formulas for the visco-elastic con- 
stants of the uniform electron liquid [23], can consider- 
ably improve the agreement between the calculated and 
the measured resistivity. 

We start by writing down the classical (single-particle) 
formula [l| for the resistivity of an electron gas of density 
no with impurities randomly distributed with density nc 



p = k F ni(Ttr{k F )/(enQ). 



(1) 



Here k F is the Fermi wave- vector of the electron gas, 
o~tr(kp) is the transport cross-section of an electron at 
the Fermi level scattered by the potential of an individ- 
ual impurity, and e is the absolute value of the electron 
charge. The basic assumptions underlying Eq. ([TJ arc: 
(i) Electrons do not interact with each other while being 
scattered by the impurities, and (ii) Electrons feel only 
one impurity at a time, i.e., the coherent scattering of 
an electron from more than one impurity is neglected. 
Both assumptions will be relaxed in the treatment that 
follows. 

Let us consider a monochromatic and uniform external 
electric field E ext (t) — E> ext e~ lut applied to electron gas 
with impurities positioned at Rfc, k = 1,2,.... We can 
write the current density averaged over the normalization 
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volume V as [? ] 

%c f 

ji{u) = —y J drdr'xij(r,r',L))E, 



extj i 



(2) 



where Xy(r, r',u>) is the current-density response func- 
tion of the inhomogeneous electron gas with impurities. 
A summation over the repeated Cartesian index j is im- 
plied. We transform Eq. ((2J with the help of the sum-rule 

r? i 



c [ui — u>„ 



4) J Xij(r,r',u)dr' = -^x 

/2 
X %k (v, r» V^V^o(r') dv' + —n Q (r) (3) 
J m 



where 



V (r) =^ Uo (r-R fe ), 



(4) 



Vq(t) is the bare potential of one impurity centered at 
origin, no(r) is the ground-state electron density, oj p — 
^/47Te 2 no/m is plasma frequency of the homogeneous 
electron gas without impurities, and c and m are the 
speed of light in vacuum and the mass of electron, re- 
spectively. Applying Eq. (j3|) twice with respect to the 
integration over r and r' in Eq. ([2]), and using the ex- 
pression for the density-response function 

Q 

X(r,r ,w) = ^Vj • Xij(r,r ,w) ■ V', 

together with the static sum-rule [l2| 

J x(r, r', 0)V^ o (r') dv 1 = V,n (r), (5) 
we eventually write the current-density as 

1 



m(uj 2 — uj 2 ) 



n E, 



ext,i ~ 



m(oj 2 — UJ 2 )V 



dr / dr' 



W v 

x[V i Vo(r)][x(r,r / ,o;)-x(r,r , ,0)][V^o(r')]]^}G) 

Independently of the foregoing considerations, just 
from Maxwell's equations for the electromagnetic field, 
we find that the resistivity can be written as [? ] 



p(uj) 



Ej(u) 
3i{u) 



— [ Anie - 
uj 



"E 2 ext 



jiEe 



(7) 



From Eq. ^ we conclude that the expression in the 
parentheses of Eq. ([7]) is zero at uj = 0. Taking the limit 
uj — » in Eq. J7| by L'Hopital's rule and using Eq. © 
again, we arrive at 
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en\V 



[W (r)-Eaf][V'7o(r')i Hi ] 



9Imx(r, r',uj) 



duj 



drdr', 



where E ex4 is the unit vector parallel to E ext . 

Equation © is the formal solution to the problem of 
expressing the resistivity in terms of the density-density 
response function x °f the interacting inhomogeneous 
electron gas with impurities. Using the relation [6| 



X 1 (r,r',u)) =x ir 1 s (r,r',o;)-/ xc (r,r / ,a;) 
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we can conveniently rewrite Eq. (J8J) in terms of the 
Kohn-Sham (KS) density-density response function Xks 
of non-interacting electrons and the dynamical exchange 
and correlation kernel f xc [? ] 



P = Pi + P2, 



(10) 



Pi 



en\y 



[VV KS (r) ■ E ext ][VV KS (r') ■ E ext ] 



dlmxKs{r,r',uj) 



duj 



drdr' , 



(11) 



where Vks(t) is the static KS potential, and 



P2 



en\y 



[V r n (r)-E ext ][V r m (r')-E e 



9Im/ xc (r, r',uj) 



duj 



drdr'. 



(12) 



The first equation (fTTjl is the single-particle (KS) con- 
tribution to the resistivity. The second (fTS]) is the dynam- 
ical exchange-correlations contribution. If the frequency 
dependence of f xc is neglected - as one does, for exam- 
ple, in the adiabatic approximation to TDDFT - then 
Eq. (HU) yields p 2 = 0. 

To establish the connection between p\ and the clas- 
sical potential-scattering result of Eq. (p}, we must ne- 
glect in Eq. (fTTj) the coherent scattering from multiple 
impurities. To do this we replace the full KS potential 
VKs{r) by the KS potential associated with a single im- 
purity in the electron gas, and we interpret the KS re- 
sponse function XKs{ r , r ',w) accordingly. The normal- 
ization volume is taken to be equal to the volume per 
impurity, i.e., V — l/n,. It can be rigourously proved [? 
] that Eq. (fTTj) . thus modified, is equivalent to Eq. (fTJ. 
This result, combined with the discussion of the previous 
paragraph, leads to the important conclusion that the 
adiabatic approximation to TDDFT is equivalent to the 
classical potential-scattering (T-matrix) approach as far 
as the calculation of the resistivity is concerned. 

The single-particle contribution to the resistivity is 
conventionally obtained from Eq. (UJ) [? ], using the T- 
matrix (phase-shift) technique to calculate the scatter- 
ing cross-section from the static KS potential @, 0] • To 
find the many-body contribution to the resistivity from 
Eq. (|T2|) . we need a good approximation to the dynami 



cal exchange and correlation kernel f xc . It is known [24 1 



(°) that f xc (r,r') is strongly non-local (i.e. a long-ranged 
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function of |r — r'|) and this non-locality is crucial to 
a proper description of many-body effects in transport 
phenomena, even on a qualitative level jl3| . This imme- 
diately poses the problem of constructing a reasonably 
accurate non-local approximation for f xc . In a recent 
paper [lj| we have shown how this can be done start- 
ing from an exact representation of the scalar f xc kernel 
in terms of the tensorial exchange and correlation kernel 
f xc of time-dependent current density functional theory. 
This representation reads 

,2 r , * r , * ,-i 



hc= V- 

c 



Xks 



7«c) f( 



. Xks' 



-fxc) T 



(xKs-hc) -Xks (tXksT) 1^} • VV 



(13) 



where xks is the KS current-density response function 
and T is the projector operator onto the subspace of 
transverse vector fields (i.e. divergence-free fields) [? 
]. By making use of the local density approximation 
(LDA) for the tensorial f xc in the right-hand side of 
Eq. (1131) . we obtain a non-local approximation for the 
scalar f xr , which satisfies the zero-force sum-rule require- 
ments [24lll3| , and can, therefore, be considered a promis- 
ingly accurate approximation for transport problems. 
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FIG. 1. (Color online) Residual resistivity of aluminum due 
to substitutional impurities of atomic number Z, as a func- 
tion of Z. The chained curve with circles (red online) is our 
result with inclusion of the dynamical exchange and correla- 
tions (the sum of pi and p2 obtained with use of Eqs. (0 and 
(112[l . respectively). The chained curve with squares (black 
online) is the result of the single-particle theory (pi only) . 
Solid squares are experimental data compiled from Refs. [251 ] . 

In Refs. [lo|, [13 , the LDA to the exchange and cor- 
relation kernel f xc of the TDCDFT has been worked 
out within the framework of the hydrodynamics of in- 
homogeneous viscous electron liquid. We, therefore, use 
Eq. (fl~2"T) with f xc given by Eq. fp~3|) and f xc as expressed 
in Ref. [26| through the viscoelastic constants of elec- 
tron liquid. In Fig. [U we present results for resistivity 



for substitutional impurities of atomic number Z from 
11 through 32 in an aluminum host. The latter is mod- 
eled as a jellium with Wigner-Seitz radius r s — 2.07. 
In this calculation we have neglected the coherent scat- 
tering from multiple impurities, focusing instead on the 
many-body dynamical exchange and correlations effects. 
The values of the viscoelastic constants were taken from 
Ref. [23| ■ Our purpose is not to take into account all the 
effects that could possibly contribute to the resistivity in 
a real solid aluminum, but rather to show that the many- 
body viscosity corrections are sizeable and indeed of the 
right order of magnitude to account for the observed dis- 
crepancy between available theoretical calculations and 
experimental data. 

The single-particle contribution p\ calculated from 
Eq. ([1} and represented by the chained curve with squares 
(black online) is found to be in agreement with earlier cal- 
culations The total resistivity, including dynamical 
exchange and correlation contributions, is represented by 
the chained curve with circles (red online) . An improved 
agreement between theory and experiment can be clearly 
seen from the figure. The effects left out by our calcu- 
lation that could possibly contribute to the remaining 
disagreement between the theory and experiment are (i) 
the band structure and lattice distortion effects, (ii) the 
possible spin-polarization, and (iii) the coherent scatter- 
ing by the impurities at different sites. Another poten- 
tially important source of error can be in the values of 
the visco-elastic constants of the electron liquid. 

Finally we note that Eq. (j8]) allows us to establish a 
general relation between the impurity resistivity and the 
friction coefficient [? ] of the same host for the same 
type of impurity atom. The latter can be written as 12 1 



Q=-J [VVb(r)-v][V'y (r')-v] 
<9Imx(r, r', lu) 



dr dr' , 



j=0 



(14) 



where v is the velocity of the atom, and, comparing with 
Eq. ([8} , we can write 



p = riiQ/(enl). 



(15) 



We point out that the relation (fl"5)) quite generally holds 
within the many-body theory and is a stronger statement 
than pi — niQi/ '(ertg), which is a simple consequence of 
Eq. ([I} and the corresponding single-particle result for 
the friction coefficient Qi = no kirat r (kp) [2t| . 

In conclusion, we have developed the non-adiabatic 
time-dependent density functional formalism for a sys- 
tematic calculation of the dc residual resistivity of met- 
als with impurities. The contribution to the resistivity 
arising from the many-body interactions has been ex- 
pressed through the dynamical exchange and correlation 
kernel f xc . We have shown that all the dynamical effects 
of the electron-electron interaction are contained in the 
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frequency dependence of f xc . The adiabatic approxima- 
tion, which neglects this frequency dependence, is exactly 
equivalent to the conventional single-particle potential- 
scattering theory of the resistivity, provided the coherent 
scattering from multiple impurities is neglected as well. 
Our calculations of the residual resistivity of Al with var- 
ious impurity atoms of different nuclear charge show that 
the inclusion of dynamical exchange and correlation con- 
siderably improves the agreement between theory and 
experiment. 
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